Higgs boson masses in an extension of the MSSM with vector-like quarks 
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We study the Higgs sector of the minimal supersymmetric standard model extended with vector- 
like quarks, at the one- loop level. The radiative corrections to the tree- level masses of the scalar 
Higgs bosons are calculated by including the contributions from the loops of top quark, vector-like 
quarks, and their scalar superpartners, for a reasonable parameter region. We find that the mass of 
the lightest scalar Higgs boson at the one-loop level should be larger than 85 GeV, if we take into 
account the negative experimental result for the Higgs search at LEP2. As the radiative corrections 
are calculated in some detail, we also find that the mass of the lightest scalar Higgs boson at the 
one-loop level is bounded from above at 280 GeV, This upper bound is increased from a previous 
result. It may provide a wider possibility for the future collider experiments to discover the lightest 
scalar Higgs boson of this model. 

PACS numbers: 



I. INTRODUCTION 

Although the minimal supersymmetric standard model 
(MSSM) [1-5] is the simplest version among the super- 
symmetric extensions of the Standard Model (SM), with 
just two Higgs doublets, there are some arguments that it 
should be extended [6] . One of the motivations to extend 
the MSSM may be the fine-tuning problem, which states 
that the radiative corrections to the tree-level mass of 
the lightest scalar Higgs boson in the MSSM should be 
large enough to be consistent with the experimental lower 
bound on the Higgs boson mass at LEP2, around 114 
GeV. If the radiative corrections are not large enough, 
the lightest scalar Higgs boson in the MSSM might face 
the difficulty of incompatibility with experiments. 

In the next-to-MSSM [7-10], where an additional Higgs 
singlet is introduced to the MSSM, the fine-tuning prob- 
lem may be alleviated by the contributions from the 
Higgs singlet. The presence of the Higgs singlet increases 
the upper bound on the tree-level mass of the lightest 
scalar Higgs boson in the next-to-MSSM up to about 
115 GeV, and thus the radiative corrections need not be 
large enough to be consistent with LEP2 data [11-38]. 

The burden of the fine tuning in the MSSM might also 
be reduced by adding extra quarks, such as the sequential 
quarks or the vector-like quarks. If the additional quarks 
are those of the fourth generation, with masses at about 
400 GeV, it has been shown that the mass of the lightest 
scalar Higgs boson in the MSSM might increase up to 400 
GeV [39]. If the additional quarks are vector- like quarks, 
as proposed about two decades ago by Moroi and Okada 
[40,41], the fine tuning problem might be solved by the 
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radiative corrections due to the contributions from the 
vector-like quarks and the vector-like scalar quarks. The 
vector-like quarks and their superpartners increase the 
mass of the lightest scalar Higgs boson well above the 
LEP2 bound, without requiring a large SUSY breaking 
scale as 2 TeV [40-46]. 

Recently, in the extension of the MSSM with vector- 
like quarks, Martin has investigated the one-loop cor- 
rections to the mass of the lightest scalar Higgs boson 
[46]. He has shown that the radiative corrections due to 
the vector-like quarks and vector-like scalar quarks may 
be as large as about 55 GeV, if the mass of the light- 
est scalar Higgs boson before correction is 110 GeV. His 
result has been obtained by neglecting some soft SUSY 
breaking parameters to simplify the 4x4 mass matrix of 
the vector-like scalar quarks into two 2x2 submatrices, 
and assuming that the weak eigenstates and the mass 
eigenstates of the vector-like quarks are identical. 

In this article, we show that the one-loop mass of the 
lightest scalar Higgs boson may be even as large as 280 
GeV if we carry out the calculations without the above 
simplifications assumed in Ref. [46], by taking into ac- 
count the contributions from top quark and vector-like 
quarks and their scalar superpartners at the one-loop 
level. The masses of the vector-like quarks are assumed 
to lie in the range of 300 GeV to 550 GeV, and the masses 
of the vector-like scalar quarks in the range of 200 GeV 
to 910 GeV. The lower bound on the mass of the light- 
est scalar Higgs boson in this model is determined such 
that it should not have been discovered at LEP2 with a 
mass smaller than its lower bound. The negative experi- 
mental result for the Higgs search at LEP2 set the lower 
bound on the mass of the lightest scalar Higgs boson as 
85 GeV. We may note that, since the mass of the light- 
est scalar Higgs boson can be sufficiently large, there are 
wide regions in the parameter space of the extensions of 
the MSSM with vector-like quarks which shall be exam- 
ined by the Higgs search in future high-energy collider 
experiments. 
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II. HIGGS BOSONS 



m\ cos P sin /3 — m'^d 



(4) 



Martin have investigated several versions of exten- 
sions of the MSSM with the vector-like matters in Ref. 
[46]. To be specific, we study the QUE model among 
them, where three vector-like chiral fields are introduced 
to the MSSM. These vector-like fields transform under 
SU{i)c X SU{2)l X C/(l) as 



Q = (3,2,1/6), 
W = (3, 1,2/3), 
f = (1,1,-1), 



Q=(3,2,-l/6), 
= (3, 1,-2/3), 
f = (1,1,1), 



(1) 



where Q, Q, U, and U are the vector-like quark fields, 
while £ and £ are the vector-like lepton fields. In terms of 
physical particles, this model has, in addition to the SM 
quarks and leptons, Vi and V2, quarks with charge +2/3, 
_B', a quark with charge —1/3, and a charged lepton L', 
together with their scalar partners Vi {i ~ 1,2,3,4), Bi 
{i = 1,2), and L, {i = 1,2). 

The superpotential of our model may be written as 

W = Wmss\i + MqQQ. + MuUU + Me££ 

+ kHuQU - hTidQU , (2) 

where Mq, M[/, and Me are the soft SUSY breaking 
masses of the vector-like sector, and k and h are Yukawa 
coupling coefficients to the Higgs fields Hu and Hd, with 
the weak hypercharge -1-1/2 and —1/2, respectively. 

The tree-level Higgs potential of our model is the same 
as the Higgs potential of the MSSM at the tree level. It 
is given as 



Vo 



mi\HuV + m'd\HdV - [midHuHd + H.c, 
1 



0(51+52') \Hu\'-\H, 



(3) 



where m^, m^, and rriud arc the mass parameters, ijj = 
{Hl,H-) and = {H+,H^^) are two Higgs doublets, 
and gi and g2 are respectively the gauge coupling coeffi- 
cients for U{1) and SU{2). There are two neutral scalar 
Higgs bosons and a neutral pseudoscalar Higgs bosons in 
the MSSM. The vacuum expectation values (VEVs) de- 
veloped by two neutral scalar Higgs bosons are denoted 
as Vd ~ (Hd) and w„ = {Hu)-, and their ratio is repesented 
by tan/3 = v^/vd- In this article, we would not consider 
the CP mixing between the scalar and pseudoscalar Higgs 
bosons. Thus, the vacuum expectation values of the neu- 
tral Higgs fields as well as the parameters of the Higgs 
potential are all real. 

Two mass parameters and can be eliminated by 
the minimum equations with respect to the scalar Higgs 
fields. Thus, just two paramaters, rriud and tan (3, de- 
scribe the tree-level masses of the scalar Higgs bosons. 
They are given by the eigenvalues of the tree-level 2x2 
mass matrix, M", for the scalar Higgs bosons. The ma- 
trix elements of M'^ are given explicitly as 



m% cos^ /? + m^d tan (3 , 
m% sin^ 13 + m^^ cot f3 , 



where = (.g^+.g|)i'^/2 is the squared mass of Z boson, 
with V = \/v'd + = 175 GeV. Note that the tree- level 
mass of the lightest scalar Higgs boson is smaller than Z 
boson mass. 

In order to study the radiative corrections to the tree- 
level Higgs boson mass, we evaluate the Higgs potential 
at the one-loop level. The one-loop effective potential is 
given by the effective pottential approximation as [47] 



niMf 
647r2 



(5) 



where A is the renormalization scale in the modified 
minimal subtraction scheme, A4i are the field-dependent 
masses of quarks and scalar quarks, and n; are the de- 
grees of freedom arising from color, charge, and spin fac- 
tors of the particles in the loops, hence, ni — —12 for 
quarks and n/ = 6 for scalar quarks. We take into ac- 
count top quark, vector-like quarks and their scalar su- 
perpartners, as their contributions are significant for the 
mass of the lightest scalar Higgs boson at the one-loop 
level in our analysis. The full Higgs potential at the one- 
loop level is thus 



V^Vo + Vi 



(6) 



The masses of the scalar top quarks after electroweak 
symmetry breaking may be expressed as 



nif r 



with 



2 / X 2 

2 I 



+ mf [ At ~ fi cot f3 



(7) 



(8) 



where rriQ and itit are soft SUSY breaking masses, 
rut = hfVu is the mass of top quark, with ht being the 
Yukawa coupling coefficient for top quark. At is the trilin- 
ear SUSY breaking mass parameter, and ^ is a parameter 
with mass dimension. Note that Xt represents the mix- 
ing between the left-handed and the right-handed scalar 
top quarks. 

Next, let us study the vector-like sector of our model. 
We adapt the notation of Re. [46]. The mass matrix for 
the vector-like quarks at the tree level is given by 



with 



MfM'f 
M^'pMf 



Mq kvu 
hvd Mu 



(9) 



(10) 



The squared masses of the two vector-like quarks, 
(i = 1,2), are given by the doubly degenerate eigenvalues 
of m%. 
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The mass matrix of the vector-hke scalar quarks is given by 

ml 6$ AkVu - k^ivd 



V \ AhVd - hy.Vu ml I ' ' 

AkVu - kfivd 



where 6$ and are the soft SUSY breaking parameters 
with mass-square dimension, mi, TO2, ''ns and TO4 are 
the soft SUSY breaking masses and and Ah are the 
trihnear mass parameters. We assume that they are all 
real and yield real eigenvalues. We note that, unlike Rcf. 
[46], Af|, cannot decomposed into two 2x2 submatri- 
ces, since we keep non-zero off-diagonal elements in m|n. 
Moreover, we do not neglect 6$ and b^. This is the main 
difference from Rcf. [46] . 

Now, we apply the above formulae for the tree-level 
masses of top quark, vector-like quarks, and their scalar 
superpartners into the one-loop effective Higgs potential 
in order to evaluate their contributions to the mass of 
the lightest scalar Higgs boson. The masses of two scalar 
Higgs bosons at the one-loop level are given as 



Tr U/ T 



Tr{M)j 



4det M 



(12) 

where Si is the lighter scalar Higgs boson as we define 
< rns2, and M is the 2x2 squared mass matrix for 
the scalar Higgs bosons at the one-loop level. 
It is convenient to decompose M as 



(13) 



where M° is the tree- level mass matrix, and Af*, , 
and respectively denotes the radiative corrections 
due to the contributions from the loops of top quark 
and scalar top quarks, from the loops of the vector- 
like quarks, and from the loops of the vector-like scalar 
quarks. Note that the factor of 2 on shows that two 
vector-like quarks produce the same radiative corrections 
due to the degeneracy between them. Let us calculate 
each of these radiative corrections. 

The radiative corrections due to the loops of top quark 
and scalar top quarks are given as = 1, 2) 



M, 



iWlW] gjml^ml) 



3A*A* 
' ^ loe 



mi TO- 
tl 12 



327r2w2 



A4 



327r2w2 



{WlA] + A\W^j) 



\og{m1Jm\) 

(to- — TO? ) 



(14) 



where 
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^11 



^22 



167r2i;2 
2 on* 



m^^At COS ip 
sin^ j3 cos P 



and 



sin^ I3D 
cos^/SD* , 
- cos /3 sin , 

A{ = 0, 

wl 



(16) 



Wl 



2m[ 
sin/3 

2to2/x(/x cot /3 - At) 

sin (3 ' 
2mlAt{At - ncot P) 
sin/3 



(17) 



with 



I [ml, ml) 



[ml - ml) 



i^log^ 



TOJ 

A2 



TO loE 

y ° fy2 



(18) 

The radiative corrections due to the loops of the vector- 
like quarks are given as (i, j = 1, 2) 

iWYWj a(m\,m\) 



327r2w2 (^2 



-iAYAX 



V2 



327r2'y2 
3 



log 



- to2,J2 

fny^ni'^2 



A4 



327r2u- 



■{wrA 



V aV 
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aYw"" 



log(TO^^/TO^J 

(to^ 



V2 



(19) 



where 



AX 



vvdh^iM^ 
vVuk'^ {v^k'^ — v^h^) 



+ 5vih' - vik') 



(20) 

where TOy. {i = 1,2) are the masses of the vector-like 
quarks. 

The radiative corrections due to the loops of the vector- 
like scalar quarks are given as — 1, 2) 



/ 2 2\ 

g[m^,m ) = 

" mi. 



7 2 
mt. TO,, , ^ 

:flog-f+2, (15) 
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A 3to| 
'J ^ 327r2 

k=l 



los 



A2 



1 



dSidSi 
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^ 327r2 
fc=i 



log 



A2 



dSi 



where the first-order derivative 9m ~ /dSi is given exphc- 
itly by 



dSi 



Aim 



Cim\, +Di 

Vk 



Vfc 



3Ami 



Vk 



C 



(22) 



and the second-order derivative d^m^y /dSidSj by 



4. 



Vk 



Vk 



D, 



dS,dS, 



^<k 



+ 



Vk 



3Ami +2Bm% +C 

Vk Vk 

B,ml^+C,ml^+D,) 



2Bn4^ 

Vk 
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■C)2 



(4to?^ 
X (3Aim|; 



3Ami +2Bml + CV 

Vk Vk 



(12m4 



2B,ml + 

Vk 

+ 6Aml^ 



2B) 



(4to^^ 4 
X {A.m^^^ 

Vk 



3Ami 

Vk 



2Bml 

Vk 



C)3 

Bim^ +Ciml +A) 



X {AjiTi 



BiTn 



Vk 



Vk 



Dj) (23) 



The exphcit expressions for the various coefficients in the 
above formulae are given in the Appendices: A, B, and C 
in Appendix I, Ai, Bi, d, and Di (i = 1,2) in Appendix 
II and Aij, Bij, Cij, and D^, ~ 1,2) in Appendix 
III. 



III. NUMERICAL ANALYSIS 

For numerical analysis, we simplify our formulae by 
assuming that k = h, Mq = Mu, mg ~ my, rni = 
= TTT-s = m-i, and 6^ = 6$. We set mi — m2 ~ 
TO3 — 7714 as M^. Note that these simplifications reduce 
the number of independent parameters, but we are still 
left with a number of them: rudu and tan/3 in M°, ^, 
toq = rriT, and At in A/*, and M^, k, and 6^ in 
and . We set the allowed ranges for these parameters 
as follows: < man (GeV) < 200, 2 < tan/3 < 30, 
150 < n (GeV) < 500, 500 < mg (GeV) < 700, 500 < 
At (GeV) < 700, < fc < 0.5, 300 < (GeV) < 
500, 300 < Ak (GeV) < 500, and 300 < (GeV) < 
500. Here, the lower bound on the value of /i is set by 
the experimental constraint on the chargino system, the 
upper bound on the Yukawa coupling coefficient of k is 
set as 0.5 in order to escape the Landau pole at a high 
energy scale. The upper bound on the SUSY breaking 
soft masses on the scalar top quark sector arc set as 700 
GeV because we do not want a large fine tuning. For the 
top quark mass, we take 171 GeV and assume that it is 
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FIG. 1: The distribution of 50,000 points of (msj, mg^), 
at the one-loop level. The allowed ranges of the parame- 
ter values are < md^ (GeV) < 200, 2 < tan/3 < 30, 
150 < /i (GeV) < 500, 500 < tuq (GeV) < 700, 500 < 
At (GeV) < 700, < fc < 0.5, 300 < (GeV) < 500, 
300 < Au (GeV) < 500, and 300 < Jh^ (GeV) < 500. 



smaller than the masses of the scalar top quarks. It is 
obtained that the masses of the vector-like quarks lies in 
the range of 300 GeV to 550 GeV, and the masses of the 
vector-like scalar quarks in the range of 200 GeV to 910 
GeV. 

Now, let us study GzzSi (* = 1)2), the coupling coef- 
ficients of the neutral scalar Higgs bosons to ZZ. In the 
SM, the corresponding quantity is given as 



GzzH = g-imz/ coaOw 



(24) 



where 9^ is the weak mixing angle, and H is the neu- 
tral scalar Higgs boson of the SM. We introduce the nor- 
malized coupling coefficients GzzSi ^ which arc defined 
as GzzSi/GzzH- These normalized coupling coefficients 
satisfy a sum rule: 



E 



^2 



1 . 



(25) 



For a given set of parameter values, we calculate the 
masses of the neutral scalar Higgs bosons at the one-loop 
level, as well as their normalized coupling coefficients. 
The LEP2 collaborations have established the model- 
independent upper bound on the squared coupling coeffi- 
cient of a neutral scalar boson to a pair of Z bosons, as a 
function of the mass of the neutral scalar boson [48] . For 
given mass of 5*1, we may compare GzzSi with the LEP2 
result. If Gzzsi is larger than the upper bound set by the 
LEP2 result, we should increase the mass Si until GzzSi 
becomes smaller than the LEP2 result. In this way, the 
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lower bound on the mass of iSi can be established. We 
examine 50,000 sets of parameter values, and select those 
sets that satisfy the LEP2 constraint. Those parameter 
sets yield the masses of the vector-like quarks within the 
range of 300 to 550 GeV, and the masses of the vector-like 
scalar quarks in the range of 200 < my (GcV) < 910 
(* = l-4). 

Our results are shown in Fig. 1, where wc plot 
{ms^Trns^) for the selected parameter sets. The lower 
bound on ms^ is determined by the LEP constraint. 
One may note in Fig. 1 that the mass of 5*1 at the 
one-loop level may be as large as 280 GeV. This upper 
bound is reasonably larger than the previous result. In 
the previous investigation, the mass of the lighter neu- 
tral scalar Higgs boson is predicted to be less than 200 
GeV [46]. Thus, the upper bound is considerably im- 
proved. We attribute this improvement to the differ- 
ence in approximations. In Ref. [46], it is assumed that 
Mq,Mu ^ kvu,hvd and = = 0. Thus, the mass 
matrix of the vector-like scalar quark is decomposed into 
two 2x2 submatrices. In our numerical analysis, these 
approximations are removed though some other simplifi- 
cations are still used. In other words, wc explore a wider 
region in the parameter space. We note that it is possible 
to examine the full parameter space of our model if no 
approximation is made in the numerical analysis. Nev- 
ertheless, ^ 280 GeV can relieve the experimental 
pressure on our model in the sense that it allows the fu- 
ture collider experiments a wider chance to discover the 
lightest scalar Higgs boson of our model. 

IV. CONCLUSIONS 

We consider the Higgs sector in an extension of the 
MSSM with the vector-like quarks and scalar quarks. A 
reasonable parameter region is set, and the neutral scalar 
Higgs boson masses are calculated using the effective 
Higgs potential at the one-loop level, where the radia- 
tive corrections due to the contributions from the loops 
of top quark, scalar top quarks, as well as the vector- 
like quarks and vector-like scalar quarks, are taken into 
account. We find that the contributions from the vector- 
like quark and scalar quark loops arc significantly large, 
if we calculate in detail, with fewer approximations, for a 
wider parameter space. The mass of the lightest neutral 
scalar Higgs boson of our model at the one-loop level lies 
within the range of 85 to 280 GeV, while being consistent 
with the LEP 2 constraint. 
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Appendix I 

The coefficients A, B and C in formulae (22) and (23) 
are expressed as in the following: 

A = A + A, 

where 

A = -2Ml^ - 2h\l 

A = A{vd O Vu,h ^ k, Ah 4-^ Ak), 

and ApQ^^ApQ + Ml 
where 

B = 3M^ + QMlMl + 6/1^2^2 ^ ^ ^^4^4 _ ^2 
-Alvl + Ah^M^vl + 2h'k\l^ - h^fi^l 
-2hkM^Vdu + 2hiiAhVdu, 

B = B{vd ^ Vu^h k,Ah ^ Ak), 

where Vdu = (u^/2)sin2/3. 

C^C + C, 

where 

C = -2M§ - 6MlM^ - 6h'vjM^ - 6M^M^ 
-2h\tMl + 2AlvlMl - 8h'My,Ml 
-Ah^k^j^M^ + 2h^n^vlM^ + AhkM^VduMl 
-ihfiAhVduMQ - 2M^ + 2h^k^fi^vj + 26^M|^ 
-2h'My^ + 2h^byd + ^AlMy^ + Ahktib^M^, 
-ihAhb^M^a + 2h^Alvj^ - 2h^eMlvj^ 
-2h^k\jvl + 2h^fi^Mlvl + AhkM^Vdu 
-Ah^iAhM^Vdu + 4:h^ nb^M^Vdu - 4,hAkb^M^Vdu 
+Ah^kMlvlvu - ^h?k^lAkvlvu, 

e = C{vd o u„, /i o fc. Ah ^ Ak). 

Appendix II 

The coefficients At, Bi, Ci, and Di (i = 1, 2) in formulae 
(22) and (23) are expressed as in the following: 

Ai = -4h^vd, 

A2 = Ai{vd ^ Vu,h ^ k,Ah ^ Ak), 
Bi = Av^h^ + 8k^Vdvlh^ + l2M^Vdh'^ 

+8M^Vdh^ - AkM^Vuh + 2nAhVuh 
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-2fcV^Wd - 2Alvd + 2kfiAkVu, 

B2 = Bi{vd ^ Vu,h <-)■ k,Ah ^ Ak), 

Ci = -SM^v^h'^ - Sk'^v^vlh'^ + UkM^vjvuh^ 
-Ak'^Vdvih'^ + SfcV^i'rf^^ + "^Alvdvlh^ 
-Wk^M^Vdvlh^ - SeMlvdvlh" + Ablvdh" 
-UM^M^^Vdh' - AMl^Mlvdh^ + Ak^Mlvlh 
-Uk/jAkvjvuh'^ + A^ih^M^Vuh^ - Ak^nAhvlh 
+8kfib^M^Vdh - SAhb^M^Vdh - 4fj,AhM^^Vuh 
+8kM^^M^Vuh - AAkb^M^Vuh + Ak^Alvdvl 
+Ak^^?Ml^Vd + ^AlMl^Vd + Ak^fib^M^Vu 
-Ak^AkM^^Vu - ikAhb^M^Vu, 

C2 = Ci{vd o Vu,h o fc, Ah o 

Di = -Qh^k^n^vl + lO/i^'feAiAfeU^v^ + Ah^M^vl 

+Ae^?Alvl ~ Sh^k^fi^M^^vj + ^^kixAhMlvl 
+8h^k^Mlvlvj - IGh^kub^M^vl + 2AlkVduVu 
+l2h?k^iAkMl^Vuvl - Gh^kn'^M^Vuvj 
-l2h?kMlMlvuvl + I2hk^ ^iAhMlvuiil 
-Gh^AhkM^Vuvl - Ghk^ti^AhVuvj - 2Ahkblvu 
-Uh^k'^lib^M^VuV^ + 12h^Akb^M^Vuvj 
^Ah^AlvduVuVd - 2k^fi^M^^Vd - 2AlM^^Vd 
+Ah^M^M^^Vd - 2k^Alvtvd + Ah^k^Mlvtvd 
+AknAhblvd - ^h^lMlvd - Sh^lM^Vd 
-Shkijb^M^^M^Vd + ShAhb^M^^M^Vd 
-QknAlAkVduVd + 2hk^iiAhvl - 2hiiAhAlvl 
+4hk''nAhM^^vl - 2hk''^?Mlvl + Ah^k^vl^vu 
-Ah^M^Aiy^ - 2k^AhkM^vl - Ah^k^lvduVu 
+4h'^knAkMlvl - 2h'^kiJ^Akvl - Ah^k^ iib^M^vl 
+4k^Ahb^M^vl + 2hnAhM^^Vu + 2fc/iAfcM^^?;„ 
-AhkM^^M^Vu - AhkblM^Vu - Uh^'k^MyduVu 
-^h^lib^Ml^M^Vu - Ak'^iib^M^^M^Vu 
+4kAhb^M^^M^Vu + AhAkb^M^^M^Vu 
+2h'^k'^fi'^VduVu + Sh'^k'^MQVduVu - 2hkfi^b1vu 
+12h^k^M^VduVu - Ak^AlMl^VduVu 
-Ah^AlMl^VduVu - 8h^k^ii^MlvduVu 
+8h^k^M^MlvduVu + 4:k^fiAhMlvd 
-8hkAhkMlvduVu + Ah^nAkM^ 

+8h?kAkb^M^Vd 
D2 = Di{vd^Vu,h^k,Ah^ Ak), 

Appendix III 

The coefficients Aij, Bij, Cij, and Dij, = 1,2) in 
formulae (22) and (23) are expressed as in tlie following: 

An = ~ih^, 

A22 = Aii{vd ^ Vu,h ^ k,Ah ^ Ak), 

A12 - A21 = 0, 

Bii = Uvlh^ + l2Mlh^ + 8Mlh'^ +8k^vlh^ 



-2k^l? - 2Al, 
B22 = Bii{vd ^ Vu,h ^ k,Ah ^ Ak), 
B12 = B21 = 8h^Vduk^ - 2hMlk + 2hnAh, 
Cii = -24Ml^vjh^ - 24k^vj^h'^ + 2AkMyuh^ 
-12M^h^ - AM^h^ - Ak^th^ + 4blh^ 
-16M|m|/i2 + 24:k^f/vdh^ + 4:Alvlh^ 
-lek^M^vlh^ - 8k^Mlvlh^ - 24:kfiAkVduh^ 
+8kiib^M^h - 8Ahb^M^h + Ak^^^M^^ 
+AAlMl^+Ak^Alvl, 
C22 = Cii{vd ^ Vu,h ^ k,Ah ^ Ak), 
C12 = C21 

= -IGk'^vjvuh^ + UkM^h^ - UkfiAkvjh'^ 
+A^ib^M^h'^ + 8Alvduh^ - l&k^Mlvduh^ 
-SkHl^uh^ - AtiAhMQ^h + AkMl^Mlh 
-~AAkb^M4,h + Ci2, 

Dii = -SOk^n^vjh'^ + 12k^vjv^h^ + 12M^vjh^ 

-12Alv^duh^ + 24k'^Ml^vj^h^ + 40fc/^AfeW>„/i^ 
-2Ak^Mlvdvlh^ - ASk^b^M^dh^ + '^^Ik^l 
+AfiAkMlvlh^ - Ukfi^M^Vduh^ - AblM^h^ 
-2AkMlMlvduh^ + 2AAkb^M^Vduh? 
+AM^M^^h^ - SblM^h' + AMlMl^Mlh^ 
-2Aeii^Miy^h^ + 2AknAhMlvlh^ 
+2e^j.\lh^ + Sk^M^vlh^ + Uk^M^h^ 
-Ak%vlh' - AAlMl^^h^ - 8k^y:'Mlvlh' 
+8eMlMlvlh' + SkAkb^M^^h^ + Ak^iAhbl 
+2Ak^iAkMl^Vduh^ - UAhkMlvduh^ 
-2Ak^fib^M^Vduh^ - 8kAhkMlvlh - 2AIM^^ 
+8kii^Ahkvlh + 8k^Ahb^M^vlh + Uk^fi'^Alvj 
~8knb^Ml^M^h - Uk^fi'^AhVduh - 2k^Alvt 
-Ak^AlMl^^ + Ak^Mlvth^ - 2fcV'M4^ 
+8Ahb4,Ml^M^h + 2Ak^ iMAhMlvduh 
+Ak^fiAhMlvl - UkfiAlAkVdu, 

D22 = Dii{vd ^ Vu, h^k, Ah o Ak), 

D12 = D21 

= 8h\lk'^ - 8Alvdvlk'^ + IQh^Mlvdvlk^ 

-ish^Myduk' + l2h^iAhMyde 

-12h^fib4,M^vjk^ + 2h^^^Vduk^ + Sh^M^Vduk^ 
+12h^M^Vduk^ - AhHlvduk^ ~ 8h^i?Mlvduk' 
+8h^MlMlvduk^ + lOh'^nAkvjk - hf?blk 
-2hM^^Mlk - 2hblMlk + Uh^fiAkM^yk 
-Gh^fi^Myk - Uh'M^Myk - GtiAlAkvjk 
-ShAhkM^uk + 8h^?AhkVduk + 2A\y.Vdu 
+16h^Akb^M^Vduk + 2hnAhM^^ - Ahkbl 
-QK^AhkMlvl + 12h^Akb^Myd - Ghfi^Ahvjk^ 
-AK^^ib^M^^M^ + AhAkb^Ml^M^ 

-8h^AlM^^Vdu + Sh^^iAkM^Vdu + D12, 
where Ahk = AhAk- 
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